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PCAC and coherent pion production by neutrinos 
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Abstract. Coherent 7r + and 7r° production in low energy neutrino reactions is discussed in the framework of the partially 



conserved axial vector current theory (PCAC). The role of lepton mass effects in suppressing the 7r + production is discussed. 
Instead of using models of pion nucleus scattering, the available data on pion carbon scattering are implemented for an analysis 
of the PCAC prediction. Our results agree well with the published upper limits for ty + production but are much below the 
W), recent MiniBooNE result for n° production. 
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PCAC AND FORWARD LEPTON THEOREM 

We discuss single pion production in coherent charged current (CC) and neutral current (NC) reactions e.g. Vn + 12 C - 
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C + 7T°. Our starting point is the general formula for neutrino scattering off a nucleus or nucleon at rest 



2 ■ da cc G 2 F cos 2 6 c „ Q 2 r 2 , n n 

Wdy = *** KE W [U aL + V aR + 2uV(Js] (1) 

already derived by Lee and Yang in 1962 JH for zero mass of the outgoing lepton. The momentum and energy transfer 
between incoming neutrino and outgoing lepton is given by q and v = E — E' , As usual Q 2 — —q 2 denotes the four- 
momentum transfer squaredQ. For Q 2 — > only the term containing the scalar cross section as survives. Here Adler's 
forward scattering theorem J2l based on PCAC predicts 
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For CC the limit Q 2 = cannot be reached. Therefore and for comparison with experiments we extrapolate to finite 
values of Q 2 by introducing a formfactor Ga = n^A./ (Q 2 + to a)- ^ n addition we include a correction (already contained 
in Adler's paper) due to the nearby pion pole in the hadronic axial vector current 031] 
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With Q^ in = mf,y/ (1 — y) the pion pole term vanishes for r?ty = 0, it is a lepton mass correction. 



1 Q 2 = q 2 ~u 2 \y = u/E; K = (W 2 - M^)/2M N ; u,v = {E + E' ± \q\)/2E. G F and B c are the Fermi coupling constant and the Cabbibo 
angle. 

2 U = V2f n o = 130.7 MeV 




FIGURE 1. 7r° production by neutrino scattering on carbon. Black histogram calculated from integrating (|7) over (t,y), red 
histogram from integrating © over y. A hadronic toy model is used (see text). 



COHERENT SCATTERING 

Coherent pion nucleus (wN) scattering is strongly peaked in forward direction distinguishing it from incoherent 
background. We therefore expect coherent single pion production by neutrinos to be well described by the PCAC 
ansatz. Like in the original Rein Sehgal (RS) paper Q4£| this approximation is assumed to hold also for the differential 
cross section 



dQ 2 dydt 



2tt 2 



l Q 



2 

min 



2Q 2 



2 • ) 

mm / 



Q 



2 

min 



(Q 2 +mi: 



da(TT+N -> TT+N) 

dt 



(6) 



and 



dQ 2 dydt At: 2 \q\ 
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This extension is by no means trivial, t is the four momentum transfer squared between the incoming virtual boson 
and the outgoing pion. Therefore t — cannot be reached. t min = f(Q 2 ) results in a very effective Q 2 cutoff for 
exponentially decreasing hadronic differential cross sections. The t-integral of e.g. (|6]l approaches Q only for Q 2 — > 0. 
Figure((T) shows as example 7r° production on carbon for E = 1 GeV. A hadronic toy model da/dt = aexp(-bt) with 
constant coefficients a = 3200mb/GcV , b = 40GcV is used. 

The RS paper |4|] evaluates the kinematical factor always at Q 2 = 0, i.e. Euv/\q\ —* (l — y)/y. At high energies the 
differences are negligible. At threshold they are very important. This is demonstrated in figure (|2^) using the hadronic 
toy model. The CC/NC ratio shown in figure $2&) approaches the limiting value of 2 cos 2 9c also only at high energies. 



THE ELASTIC PION NUCLEUS CROSS SECTION 



A simple model for elastic pion nucleus scattering, which can be easily implemented into MC generators is also 
contained in 01 . We discuss it here for isoscalar targets of atomic mass A leading apart from electromagnetic 
corrections to identical cross sections for 7r ±!0 . Starting from 

do(7xN — ► tcN) ^2 da c \ 
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FIGURE 2. a) Ratio cr^jj / ctJ lo of the integrated cross sections of (T7J versus the energy of the incoming neutrino. o"si mp i e is 

calculated using the kinematical approximations of 0]. b) Ratio /a" of the integrated cross sections versus the energy of the 
incoming neutrino. 



the elastic differential pion nucleon cross section at t — is calculated with the help of the optical theorem 
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where the total pion proton cross sections are taken from data. The slope &rs 1S determined via the optical model 
relation 

bns = Ir 2 A 2 / 3 e.g. R = 1.057fm . (10) 
Finally using a simple geometrical picture the absorption factor 

9^/3 

is calculated from data for inelastic pion proton scattering via 

a- , + a- , 

„ met met /n\ 
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Although this model has its limitations (e.g. it predicts da/dt—>0 for A — > oo) it has been very successful in describing 
high energy coherent neutrino scattering y| . 

In order to obtain a more precise prediction for the pion nucleus cross section in the resonance region the parameter- 
ization of the pion nucleon cross sections used in RS [4] has been replaced by detailed fits to the 7r ± p data published 
by the Particle Data Group ]5|]. An example is shown in figure (|3^). The curve labelled RS2009 in figure ([3J3) shows 
the resulting total CC coherent neutrino cross section for energies up to 2 GeV. There exist various implementations of 
the RS-model which, however, obtain different results. An example is displayed in figure (|3p). The predictions of other 
Monte Carlo generators claiming to use the RS-model show even more pronounced discrepancies [7]. The reason for 
these differences remains a puzzle. 

For the resonance region with its rapidly varying cross sections and angular distributions the hadronic RS-model is 
probably too simple. It describes badly the low energy experimental data on elastic ir 12 C scattering. Instead of refining 
it we - in the spirit of Adler's theorem - directly revert to the measured 7r 12 C cross sections. Pion carbon scattering 
data with 30 < < 776 MeV of various experiments have been subjected to a phase shift analysis and extrapolated 
to Tjr = 870 MeV by the Karlsruhe group [8]. Figure (0^) shows an example for the differential cross section da/dt 
reconstructed from these phaseshifts at = 162 MeV, close to the maximum of the first resonance. The forward 
scattering containing the bulk of the cross section is then fitted by a aexp (— bt) ansatz resulting in energy dependent 
coefficients a, b ||9|]. Using this parameterization the pion carbon elastic cross section in the resonance region is below 
the hadronic RS-model but approaches it quickly at higher p^, see figure d4j?). 
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FIGURE 3. a) Fit to the total 7r + p cross section using the PDG tables [5] for pion laboratory momenta up to 2 GeV. b) Total CC 
cross section versus neutrino energy using the updated hadronic RS model (RS2009). For comparison the prediction used by the 
SciBooNE Collaboration 16] is also shown. 
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FIGURE 4. a) Differential cross section da/dt for elastic pion carbon scattering. Th red line labelled 'phasehifts' is calculated 
from the phasehifts in [8]. The green line (BS2009) is an exponential fit to the forward cross section. The blue line (RS2009) 
represents the result of the updated RS hadronic model, b) Total elastic pion carbon cross section versus the laboratory pion 
momentum in the updated RS hadronic model (RS2009) and from exponential fits to the phaseshift analysis (BS2009). 



RESULTS 

Using the fits discussed in the preceding section we get a substantial modification of the PCAC prediction for pion 
production off carbon nuclei for NC and CC reactions at low neutrino energies. This is demonstrated in figure (O 
where the new results are compared with calculations using the updated hadronic RS-model. 

Our predictions for the total cross section are compatible with other PCAC based calculations [ 10] and remarkably 
close to certain variants of microscopic nuclear physics models flTl [l2tl . Differential distributions are more sensitive 
to model details. The MiniBooNE collaboration has proposed to use E n (l — cosO„.) as variable for the analysis of 
neutrino scatterin g IU3I1 (E v and 0„- defined in the laboratory system). As can be seen in in figure (J5J a recent nuclear 
physics model B14l[l2ll agrees well with our PCAC model at a neutrino energy typical for the MiniBooNE experiment. 

The extension of the new ansatz to other nuclei is of particular importance. At this moment we propose to use in the 
spirit of the optical model an A 2 / 3, scaling law which is close to the effective A-dependence obtained in the hadronic 
RS-model for light nuclei. 

Our results agree with the published experimental limits on coherent 7r + production lfl5l [l6ll . Using the A 2 / 3 
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FIGURE 5. Cross section per nucleus of coherent tt production by neutrinos off carbon nuclei, a) NC reaction + C — > 
+ 12 C + 7r°, b) CC reaction + 12 C — > fj,~ + 12 C + 7r + . The data in units of 10 -40 cm 2 are plotted versus the neutrino energy 
in GeV. The upper curve is calculated using the hadronic RS model, the lower curve using our parametrization of pion carbon 
scattering data. 
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FIGURE 6. Differential cross section da/dE^il — cosGjr) versus E v (l — cosGtt) calculated fur NC pion production at a 
neutrino energy of 0.8 GeV. The histogram uses the BS2009-model |9] and the curve represents the results of a recent nuclear 
physics calculation flUl . 



scaling law it also agrees with the n° data of the Aachen Padova experiment fl7il . Like all other recent theoretical 
predictions we have a problem with the MiniBooNE ir° result lfl3ll . PCAC models have very little flexibility in tuning 
the predictions. Before, however, claiming that the model is falsified one would like to clarify several questions, e.g. 
the dependence of the experimental result on the use of inappropriate Monte Carlo models and how the experiments 
ensure the coherence of the process. 

Instead of using a two parameter fit of the pion carbon differential cross section a variant of the new PCAC model 
has been studied in which the full angular distribution as represented by the phase shifts (see figure (0^)) is utilized. 
For energies between two data sets a linear extrapolation of the differential cross sections da/dt at a fixed scattering 
angle in the CMS system is applied. The resulting total neutrino cross section and the pion angular distribution in 
the forward direction changes only at the level of a few percent. In contrast to the simpler model there is, however, a 
long tail at larger angles. An example is shown in figure (0. These differences might become important when precise 
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FIGURE 7. Tail of the differential cross section da /dE n ( 1 — cos 6^) versus ( 1 — cos Q-k ) calculated for NC pion production 
at a neutrino energy of 1.0 GeV. The black histogram is obtained using the two parameter fit to pion carbon scattering described 
in 0]. The red histogram uses the full angular distribution. For comparison a prediction using the updated hadronic RS-model is 
also shown (green histogram). 



experimental data are available. With only a few sets of phaseshifts for other nuclei on-hand the extension of this 
model version to non carbonic targets requires further research. 
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